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TECHNICAL PUBLICATION 


USING PLATE FINITE ELEMENTS FOR MODELING FILLETS IN DESIGN, 
OPTIMIZATION, AND DYNAMIC ANALYSIS 

1. INTRODUCTION 


Fillets are one of the most commonly found design features in structures. They are used for reasons 
varying from reducing stress concentrations to facilitating machining. In many cases, a finite element 
model of the structure with these fillets is required for design and analysis. For component-level modeling, 
either plate or solid finite elements will usually be required to give the level of fidelity desired. Upon 
immediate inspection, it would appear that solid modeling would be required to model the fillets accu- 
rately. Frequently, however, the most accurate and economical element for modeling the structure is the 
plate element. The reasons for this are that plates can more accurately predict the behavior of “plate-like” 
structures (defined as the ratio of the thickness to the other dimensions being <10 percent), and that the 
reduction in the number of degrees of freedom (DOFs) in the model is usually in excess of 90 percent. Even 
considering the enormous power of present day computing packages, this reduction in DOFs over solids is 
desirable for a variety of reasons. First, the design may require many iterations of a model or numerical 
optimization. Second, model creation is significantly simplified by the use of plates compared with accu- 
rately meshing fillets using solid elements; many automeshers generate unacceptable meshes of these 
shapes, and manual meshing is nontrivial. In addition, a smaller model significantly facilitates postprocessing. 
Finally, there are many situations where memory and/or speed of the computing environment is and will 
continue to be limited, such as in educational or small business settings. 

The simplest method of using plates in these situations is simply to ignore the fillet. Obviously, the 
fillet stresses would not be obtained using this approach, but many types of analysis do not require a 
solution for detailed stresses. These include dynamic analyses, design for deflection or vibration control, 
or preliminary design. Clearly, though, ignoring fillets would cause an underprediction of the stiffness of 
the structure, especially for thin-walled configurations. Therefore, this problem provides the motivation for 
this research, which is to generate a method for using plate elements themselves to represent the stiffness of 
the fillets. This method will allow the user to take advantage of modeling with plate elements compared 
with the substantially more numerically intensive solid elements and to still accurately predict structural 
stiffness. 

One example of the type of structure this method would be applicable to is the structural platform 
for use in a fine-motion stage 1 as shown in figure 1 . The complexity of such structures leads to the obvious 
use of finite element methods for evaluating the performance of the design. Since a computer-aided design- 
based solid geometry model was created initially, the generation of a finite element mesh using solid 
elements was attempted using automatic meshing software. However, the intersection of thin-walled plates 
is problematic for these automeshers; a substantial number of tetrahedral elements with unacceptable 
aspect ratios are created at the fillets. In addition, the multifaceted intersections require an extremely large 
number of elements; design changes and subsequent analysis therefore take a considerable amount of time 
and effort. 



Figure 1. Four-inch cube, thin-walled structure 
used for 6 DOF fine motion stage, 
with wall thickness of 0.040 in 
throughout. 


Based on the geometry of these fillets, various analysis methods were examined to try to obtain an 
accurate, yet timely, modeling methodology. Plate theory was initially investigated for the analytical basis 
of this study. Liu and Chang 2 and Ganesan and Nagaraja Rao 4 have performed static and vibration analysis 
for plates whose thickness varies. Rayleigh 4 and Timoshenko 5 have also covered vibrations of plates with 
uniform thickness with closed-form solutions stated for simple cases. The most extensive collection of 
closed-form solutions of uniform and nonuniform thickness plate theory comes from Leissa. 6 Cases of 
highly complex geometries are assessed experimentally to gain insight into the behavior of a complex plate 
constrained by simple boundary conditions. Correction factors are found for the tested geometries where 
strengthening ribs or webs are used and subsequently inserted back into the basic equations for simplified 
cases. The conclusion from this study was that the closed-form solutions for this type of analysis are 
limited to very simple geometries where the boundary conditions and geometry are easily defined, and that 
for the general fillet case, the boundary conditions prohibit strict analytic solutions. 
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2. THEORETICAL DEVELOPMENT 


Several methods were assessed for achieving the goal of this research. The first proposed solution 
was to simply create a dense mesh of the plates at the filleted intersection and increase the thickness of the 
plates so they would match the basic shape of the fillet, as seen in figure 2. The mesh density is based on the 
geometry of the fillet only, and the material properties of the “fillet” plates were not altered. One weakness 
of this method is that as the length of the plate element decreases to allow accurate matching with the 
thickness of the fillet, the aspect ratio becomes unacceptable unless the depth of each element is also 
reduced, which would dramatically increase the mesh density and complexity of the entire model. In addi- 
tion, this section of increased density and thickness plates would have to be reproduced along both in-plane 
axes of the fillet. Because of these weaknesses, this method was deemed to be undesirable. 



Figure 2. Method of modeling fillets by stepping 
up thickness to replicate fillet geometry. 


After substantial experimentation with this and various other methods, some of which are docu- 
mented in previous work by the authors, 1 the procedure finally converged upon was to match the rotation at 
the tangent of the fillet to that of a “bridge” system. The rotation value was seen during the studies to have 
the most effect on accurate deflections throughout the structure. The fillet would be replaced by a system of 
two or three plate groups (see fig. 3), one set colinear with the nonfilleted section and a group in place of 
each “half' fillet section with endpoints at the fillet tangents. “Pseudo” thickness and stiffness (Young’s 
modulus (£)) properties for the bridge groups are calculated using the derivation shown below, while the 
colinear section uses the same plate properties as the nonfilleted sections. This geometrical configuration is 
relatively simple to implement and can be used for both doubly filleted and singly filleted sections. It is, 
however, limited to cases where the filleted intersection forms a 90° angle. The calculated properties were 
seen to be a function of two parameters: the ratio of the fillet radius (r) to the thickness of the plate (f) and 
the ratio of the thickness of the perpendicular nonfilleted sections to the colinear nonfilleted section thick- 
ness (t wa )|). In previous work by the authors, only the thickness of the bridge element (t b ) was generated, 
and this value was obtained by empirically matching the transverse stiffness of the bridge with a solid finite 
element model. Upon further examination and discussion, though, it was seen that the rotational stiffness 
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was left out of this derivation. Implementing this stiffness was achieved by matching not just a single 
stiffness constant but by using the 2x2 stiffness matrix of the fillet system at the tangent node to equate the 
rotations at the tangent. This procedure is outlined below. 



Figure 3. Representation of fillet with wide beams. 


First, upon further examination of a generic fillet configuration, it is recognized that the bending 
stiffness of the fillet itself is only a function of its length (L), and not of its width ( b ). This is essentially the 
definition of a structure in a state of plane strain, which allows the stiffness to be reduced to in-plane 
transverse and rotational stiffness only. Using this assumption eliminates the out-of-plane DOFs and 
allows the use of “wide-beam” theory for the stiffness matrix calculations rather than plate theory, 7 which 
greatly simplifies the problem. This theory states that for plane-strain conditions, beam-bending equations 
can be used for the calculation of deflection and in-plane rotation of plates by replacing Young’s modulus 
(E) with an adjusted value ^wjde-beanv w hich is obtained using the equation 

F - - (1) 

c ' wide-beam — i ■ 

1- V 

These theories were tested and verified for simple beams and plates using the finite element method, 
as described in section 3. 

In addition, if small angles are assumed, the only contribution to the stiffness from the bridge 
elements is axial extension, while the only contribution from the transverse element is bending. Referring 
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to figure 2, since the only DOFs of interest are the node 4 transverse deflection (v 4 ) and rotation (0 4 ), the 
stiffness partitions associated with node 4 are formulated from elements 1, 2, and 3. Using the partition for 
a rotated plane frame element from Cook, 8 the following is obtained: 
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For 0* = -45° and 0 3 =45°, elements 1 and 3 will yield identical results with Ly=T-L. It should also 
be noted that L is the length of the frame element, which will extend from the tangent point of the fillet to 
the midplane of the wall side, yielding Z, brid „ e = r+r wa n/2. If these are called “bridge” elements, and the 
partitions superimposed, the following is therefore obtained: 
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Applying the plane-strain assumption allows the use of a unit depth, which reduces the stiffness partition 
further to 
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At this point, the rotation (0 4 ) of the actual fillet to any arbitrary load (P) and moment (M) at the tangent 
point of the fillet is determined. First, for a fillet of radius (r), the height (/y) is related to the location (x) by 


hj - r- x(2r- x) , (6) 

where x is measured from the root of the fillet. The total thickness of the fillet (y) is therefore 



and the second moment of area of the fillet (If) is defined by 


(7) 
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where the base (b) is equal to the unit depth 1. The rotation 9f(x) can therefore be expressed as 


J(M + Pxr) A + c (9) 

EI f 

Since the boundary condition on the root side of the fillet is zero rotation; i.e., v'(jc=0)=0, the constant of 
integration C will equal the negative of the value obtained from the indefinite integral. The symbolic 
integration is performed using Mathematica™ 4.1 but is not included in this Technical Publication. d 4 is 
determined by setting x=r in equation (9). 

Similar expressions for the rotation at node 4 are now obtained for the bridge system. These are 
obtained by multiplying the inverse of the stiffness partition [K\ 44 by a column vector of the arbitrary 
moment and load. A solution is being sought such that the rotation is the same for the bridge system as the 
actual fillet, so the solutions are equated for this value from each system. Since the solutions have to be 
valid for any value of M or P, the coefficients for each of these parameters can therefore be collected and 
equated, leaving two nonlinear equations for the unknowns E h and t b : 
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Using Mathematica 4.1, this can then be solved symbolically for the desired parameters E h and t b in terms 
of r and r wal] . The solution was verified by substituting the values back into the bridge and fillet system 
rotation equations to check for equivalence. The solution is readily obtained using Mathematica, but for 
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ease of use, results for the most applicable range of parameters have been tabularized (tables 1 and 2) as 
well as fit to an explicit equation. The tables can be imported as matrices into the program Matlab™, where 
the “surfht(twall, r, eb)” command will yield an interactive interpolated grid where the E^E or t h /t value 
can be easily obtained. The surface fit of these data is shown in equations (11) and (12) where x=r/t and 
y =t wal \/t. The maximum errors of the equations from the data are 2.8 percent, and are generally <1 percent: 
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Table 1 . E b /E normalized bridge Young’s modulus. 
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Table 2. t b !t normalized bridge thickness. 
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Finally, for dynamic applications, it is also important to accurately represent the mass of the origi- 
nal fillet. This is accomplished by deriving a pseudo mass (p b ) that can be easily assigned to the bridge 
elements. This equation is obtained as a function of the actual density of the fillet (p^), assuming the use of 
a homogenous material throughout the fillet, and by matching the volume of the bridge elements with the 
volume of the original fillet: 
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3. VALIDATION 


Several levels of validation are required for verification of the technique described. The first is to 
verify the plane-strain assumption of a transversely loaded fillet. This assumption was verified by building 
a high-fidelity solid model of a cantilevered 1-in- (2.54-cm-) wide filleted section (radius of 0.25 in 
(0.635 cm), thickness of 0.04 in (0.1016 cm)), loading it transversely with a 1-lb (4.448-N) static load, and 
comparing the directional strains at a node at the tangent of the fillet. The results are in-plane strains of 
501 //strain and 243 //strain, and an out-of-plane strain of 1 //strain, verifying the assumption. The capabil- 
ity of wide-beam theory to represent plane-strain configurations was then examined. For a cantilevered 
plate with a 10-in (25.4-cm) length, 100-in (254-cm) width, 1-in (2.54-cm) thickness, with a 1-lb/in (1.75- 
N/cm) end load, the theoretical value of displacement (5) is 3.8037x10^ in (9.6614X10 -4 cm) using the 
“wide-beam” adjusted Young’s modulus where 


s= Pll + 6(PLl < 14 ) 

3 El 5(AG) 

A high-fidelity plane-strain MSC/NASTRAN finite element model yielded the exact same result. 

The next validation step is to verify the validity of using beam theory to generate the fillet stiffness. 
For the following parameters, 

£ = 10.6 xl0 6 psi (73.084 GPa), r = 0.25 in (0.635 cm), t = 0.064 in (1.626 cm), v = 0.33, 

P = l]k (6.894 kPa), Af = 1.867 lb-in (21.093 Nt-cm), and f wall =0.064 in (1.626 cm) , (15) 

the rotation at the tangent point was calculated using the fillet wide-beam theory outlined above to be 
5.3134X10 -4 rd. For comparison, a high-fidelity plane-strain finite element model (fig. 4) using the same 
parameters resulted in a rotation at the tangent of 5.76156X10 -4 rd, an error of 7.78 percent. This error is 
larger than anticipated, but is acceptable for use in the methodology developed here. 

Finally, a comprehensive test was performed by applying the technique on a representative struc- 
ture that could be modeled with a high-fidelity solid mesh and that could be constructed and tested. 
A 5x2.5x1.26 in (12.7x6.35x3.2 cm) double “box” structure was machined and modeled to meet these 
criteria. A photo of the structure is shown in figure 5, and the highly dense solid mesh and plate mesh using 
the bridge system are shown in figures 6 and 7. To assess the impact of the fillet, the structure was also 
modeled using plates but with the fillets completely ignored. 
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Figure 4. Fillet high-fidelity plane-strain finite element model 



Figure 5. Typical filleted thin-walled structure. 
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Comparisons of the size of the models and results of a modal analysis and test are shown in table 3. 
Mode 8 for the two models is shown in figures 8 and 9. The solid model has 5 1,520 nodes (257,600 DOFs) 
and 45,120 elements, while the plate “bridge” element model has 4,465 nodes (22,325 DOFs) and 4,392 
plate elements. The results show that the application of the technique results in at least a 90-percent reduc- 
tion in the number of DOFs with no loss of accuracy, and that ignoring the fillets underpredicts the natural 
frequencies substantially. 


Table 3. Comparison of models and test for double box. 


Free-Free 

Mode 

Number 

Test 

Frequency 

(Baseline) 

Solid 

Model 

Frequency 

Error 

(%) 

Plate Element 
Model With 
“Bridge” Fillets 

Error 

i%) 

Plate Element 
Model 

Without Fillets 

Error 

(%) 

7 

484.0 

487.6 

0.8 

494.6 

2.2 

449.1 

-7.2 

8 

540.0 

534.3 

-1.1 

541.0 

0.2 

467.6 

-13.4 

9 

632.0 

627.4 

-0.7 

633.8 

0.3 

547.1 

-12.8 

10 

1,012.0 

999.7 

-1.2 

998.7 

-1.3 

1,000.7 

-1.1 

11 

1,412.0 

1,403.4 

-0.6 

1,408.8 

-0.2 

1,282.2 

-9.2 

12 

1,900.0 

1,903.6 

0.2 

1,921.8 

1.1 

1,687.7 

-11.2 

13 

2,320.0 

2,290.8 

-1.3 

2,322.4 

0.1 

1,992.3 

-14.1 

14 

2,810.0 

2,792.6 

-0.6 

2,806.6 

-0.1 

2,291.4 

-18.5 



Figure 8. Solid model, free-free mode 8 at 524.5 Hz. 
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4. CONCLUSIONS AND FUTURE WORK 


A methodology has been developed that allows the use of plate elements instead of numerically 
inefficient solid elements for modeling structures with 90° fillets. The technique uses plate “bridges” with 
pseudo Young’s modulus ( E b ) and thickness (t h ) values placed between the tangent points of the fillets. 
These parameters are obtained by solving two nonlinear simultaneous equations in terms of the indepen- 
dent variables r/t and r wal j/t. These equations are generated by equating the rotation at the tangent point of 
a bridge system with that of a fillet, where both rotations are derived using beam theory. Accurate surface 
fits of the solutions are also presented to provide the user with closed-form equations for the parameters. 
The methodology was verified on the subcomponent level and with a representative filleted structure, 
where the technique yielded a plate model exhibiting a level of accuracy better than or equal to a high- 
fidelity solid model and with a 90-percent reduction in the number of DOFs. The application of this method 
for parametric design studies, optimization, and dynamic analysis should prove extremely beneficial for 
the finite element practitioner. Although the method does not attempt to produce accurate stresses in the 
filleted region, it can also be used to obtain stresses elsewhere in the structure for preliminary analysis. A 
future avenue of study is to extend the theory developed here to other fillet geometries, including fillet 
angles other than 90° and multifaceted intersections. 
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